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Abstract 

Binary systems containing at least one radiopulsar are excellent laboratories to test 
several aspects of fundamental physics like matter properties in conditions of extreme 
density and theories of gravitation like the Einstein's General Theory of Gravitation 
(GTR) along with modifications/extensions of it. In this Chapter we focus on the 
perspectives on measuring the moment of inertia of the double pulsar, its usefulness in 
testing some modified models of gravity, and the possibility of using the mean anomaly 
as a further post-Keplerian orbital parameter to probe GTR. 



PACS numbers: 04.80.-y, 04.80.Cc, 95.30.Sf, 97.60.Gb, 97.60.Jd 



1. Introduction 

In this Chapter we discuss some applications of astrophysical binary systems containing at 
least one radiopulsai^ to fundamental physics. In particular, in Section |2?] we investigate the 
persp ectives on measuring the pulsai^'s moment of inertia dloriol. l2009al : iKramer and Wex , 
20091) through the extension of the general relativistic Lense-Thirring effect to the double 
pulsar system J0737-3039. Section ID sh ows how the double pulsar can be used to put on 
the test some models of modified gravity (Iloriol.l2009bh . In Section|4]]we deal with the pos- 
sibility of using the mean anomaly as a further post-Kepleria n parameter u seful to test the 
gravitolectric part of the General Theory of Relativity (GTR) ( loriol . 2007a); for a proposal 
concerning a new post-Keplerian parameter t o test different aspects of gra vitomagnetism 
with respect to the Lense-Thirring effect, see (IRuggiero and Tartaglial.l2005h . 
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2. The Moment of Inertia of Pulsar J0737-3039A: Prospects for 
its Measurement 



The measurement of the moment of inertia / of a neutron star at a 10% level of accuracy 



(EOS) describing matter inside neutron stars ( 


Morrison et al. 


Lattimer and Schutz. 


2005; 


Lavagetto et al. , 


2007). 



After the discovery of the double pulsar PSR J0737-3039A/B system (IBurgay et al. 



2003h . whose relevant orbital parameters are listed in Table[Il it was often argued that such 
a measurement fo r the A pulsar via the post-Newtonian gravitomagnetic spin-orbit peri - 



astron precession (|Barker and O'Conneli Il975al : iDamour and Schaefer , 



would be possible after some years of accurate and continuous timing. 



19881: Iwexl 



Lvne et all (l2004l) 



1995 ) 



Table 1. Relevant Kepl erian and pos t-Keplerian parameters of the binary sys- 
tem PSR J0737-3039A/B dKramer et all 1^006). The orbital period Pb is measured 



with a precision of 4 x 10~^ s. The projected semimajor axis is defined as x = 
(«bc/c) sini, where abc is the barycentric semimajor axis (the relative semimajor axis 
a = (xA + xb)c/ sin i), c is the speed of light and i is the angle between the plane of 
the sky, perpendicular to the line-of-sight, and the orbital plane. The eccentricity is 
e. The best determined post-Keplerian parameter is, to date, the periastron rate uj 
of the component A. The phenomenologically determined post-Keplerian parameter 
s, related to the general relativistic Shapiro time del ay, is equal to sini; w e have con- 



servatively quoted the largest error in s reported in dKramer et al.. 20061) . The other 



post-Keplerian parameter related to the Shapiro delay, which is used in the text, is r. 



fb (d) (s) XB (s) e a;(degyr s r(ps) 

0.10225156248(5) 1.415032(1) 1.5161(16) 0.0877775(9) 16.89947(68) 0.99974(39) 6.21(33) 



write: "Deviations from the value predicted by general relativity may be caused by contri- 
butions from spin-orbit coupling (IBarker and Q'Connelll . Il975bh . which is about an order 
of magnitude larger than for PSR B1913-I-16. This pote ntially will allow us to measure the 
moment of inertia of a neutron star for the first time ( Damour and Schaefer . ^98 8. : ,Wexi 



19951) ." 



According to iLattimer and Schutzl (|2005l ). "measurement of the spin-orbit perihelion 
advanc e seems possible." 

In ( Kramer et al. . 20061) we find: "A future determination of the system geometry and 
the measurement of two other PK parameters at a level of precision similar to that for 
u, would allow us to r neasu r e the mome nt of inertia of a neutron star for the first time 
(IDamour and Schaefen. Il988l : IWexl 119951) . [...] this measurement is potentially very dif- 
ficult [...] The double pulsar [...] would also give insight into the nature of super-dense 
matter.'' 

In (IDamoun. l2007h it is written: "It was then pointed out (IDamour and Schaefer!. 1 1 98 8h 
that this gives, in principle, and indirect way of measuring the moment of inertia of neu- 
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tron stars [...]. However, this can be done only if one measures, beside^ k, two other PK 
parameters with 10^^ accuracy. A rat her tall order which w ill be a challenge to meet." 

Some more details are released by iKramer et al.l (120051) : "[...] a potential measurement 
of thi s effect allows the mo r nent o f inertia of a neutron star to be determined for the first 
time (IDamour and Schaefen. 1 19881) . If two parameters, e.g. the Shapiro parameter s and 
the mass ratio TZ, can be measured sufficiently accurate, an expected ujexp can be computed 
from the intersection point." 

Here we will examine with some more details the conditions which would make fea- 
sible to measure Ia at 10% or be t ter in the PSR J0737-3039A/B system in view of the 
latest timing results ( Kramer et al. . 2006h . In particular, we will show how important the 
impact of the mismodelling in the known precessional effects affecting the periastron rate 
of PSR J0737-3039A/B is if other effects must be extracted from such a post-Keplerian pa- 
rameter. Such an analysis will turn out to be useful also for purposes other than measuring 
gravitomagnetism like, e.g., putting more severe and realistic constraints on the parameters 
entering various models of modified gravity (See Section[3l. I ndeed, in doing that for, e.g., 
a uniform cosmological constant A Ijetzer and Serenol (120061) took into account only the 
least-square covariance sigma of the estimated periastron rate (6.8 x 10~^ deg yr~^): the 
systematic bias due to the first post-Newtonian (IPN) periastron precession was neglected. 
Concerning the use of the mean anomaly ^A for testing IPN effects on it in pulsar systems, 
see Section l4l 



2.1. The Systematic Uncertainty in the IPN Periastron Precession 



By assuming / w 10^^ kg m^ ( Morrison et al. . 2004; Bejger et"aD . l2005l ). the gravitomag- 
netic spin-orbit periastron precession is about wqm ~ 10~^ deg yr~^, while the error 
(^Wmeas with which the periastron rate is phenomenolq gically estimated from timing data 



is currently 6.8 x 10~^ deg yr~ (IKramer et al.l |2006|) . In order to measure the gravito- 
magnetic effect— and, in principle, any other dynamical feature affecting the periastron— 
^Wmcas is certainly of primary importance, but it is not the only source of eiTor to be care- 
fully considered: indeed, there are other terms contributi ng to the periastron precession 
(first and second post-New t onian, qu adrupole, spin-spin (IBarker and Q'Connelll . Il975al : 
Damour and Schaefeiill988l : IWexl.ll99 5)) which must be subtracted from comeas, thus intro- 
ducing a further systematic uncertainty due to the propagation of the eiTors in the system's 
parameters entering their analytical expr essions. A preliminary analysis of such aspects, can 
be found in (ILattimer and Schuta . 120051) . However, apart from the fact that its authors make 
use of the value for i measured with the scintillation method (|Coles et al. 1. 1200.5, ) which is 
highly uncertain for the reasons given below, in using the third Kepler law to determine the 
sum of the masses they also confound the relative projected semimajor axi^ a sin i (see eq. 
^) with the barycentric projected semimajor axis x, which is the true measurable quantity 
from timing data, so that their analysis cannot be considered reliable. The semimajor axis 
a of the relative motion of A with respect to B in a binary system is just the sum of the 
semimajor axes abc of A and B with respect to the system's barycenter, i.e. a = a^^ + a^^. 



*The parameter k is directly related to tlie periastr on precession lj. 

^Our a must not be confused with the one used in jLattimer and Schut3 . l2005t) which is, in fact, Obc- 
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Let us, now, consider th e largest contribution to the perias tron rate, i.e. the IPN preces- 
sion (IDamour and Deruellel 119861 : iDamour and Tayloii 1 19921) 



WlPN 



(1 



-5/3 



27r 



2/3 



(1) 



where Tq = GTXq/c^ and M = mx + m-Q, in units of Solar mass 50^0 . It is often referred 
to as gravitoelectric in the weak-field and slow-motion approximation: in the context of 
the Solar System it is the well known Einstein Mercury's perihelion precession of about 43 
arcsec cy^^. Thus, 

WgM = l^meas — ^IPN — W2PN, 

(2) 

The sum of the masses M enters eq. ([Hi; as we will see, this implies that the relative 
semimajor axis a is required as well. For consistency reasons, the values of such parameters 
used to calculate eq. ^ should have been obtained independently of the periastron rate 
itself. We will show that, in the case of PSR J0737-3039A/B, it is possible. 
Let us start from the relative semimajor axis 



a = (1 + 7^) 



CXA 



smz 



^.78949386 x 10* m. 



(3) 



It is built in terms of TZ, the projected semimajor axis xa and sin i; the phenomenologically 
estimated post-Keplerian parameter s determining the shape of the logarithmic Shapiro 
time delay can be identified with sin i in GTR and the ratio TZ = xb/xa has been obtained 
from the phenomenologically determined projected semimaj or axes, being equal to the ra- 



tio of the masses in any Lorentz-invariant theory of gr avity (|Damour and Deruellel. 11985 
Damour and Schaefen. 1 19881 : lOamour and Tayloii Il992r) 



7^ = ^ + o^4 



mB \c 

The uncertainty in a can be conservatively evaluated as 



5a < 5a\'ji + 5a\s + 5a\ 



(4) 



(5) 



with 



' Sa\n < (^) Sn = 466, 758 m, 



6a\s < a 



342, 879 m, 



(6) 



5a\xj^ < a ^ 



^ 1 = 621 m. 

XA 



Thus, 



5a < 810,259 m. 
eq. ^ yields a relative uncertainty of 



(V) 



(8) 
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It is important to note that xb, via TZ, and s have a major impact on the overall uncertainty 
in a; our estimat e has to be considere d as conservative because we adopted for 6s the largest 
value quoted in (IKramer et al.U2006h . In regard to the inclinati on, we did not use the more 
precise value for i obtained from scintillation measurements irFI( Coles et al. , 2005 ) because 
it is inconsistent with that derived from timing measurements (IKramer et all. 120061) . More- 
over, the scintillation method is model-dependent and it is not only based on a number of 
assumptions about the interstellar medium, but it is also much more easily affected by var- 
ious other effects. However, we will see that also xa has a non-negligible impact. Finally, 
let us note that we purposely linearly summed up the individual so urces of errors in vie w of 
the existing correlations among the various estimated parameters ( Kramer et al. . 20061) . 

Let us, now, determine the sum of the masses: recall that it must not come from the 
periastron rate itself. One possibility is to use the phenomenologically determined orbital 
period Pb and the third Kepler law getting 



Gm = a'^ 



2tt\- 



(9) 



With eq. ^ and eq. (O we can, now, consistently calculate eq. ([T]) getting 

2 /n_\ 3 



WlPN 



(1 



XA + Xb 



^ = 16.90410 degyr"^; (10) 



in this way the IPN periastron precession is written in terms of the four Keplerian parame- 
ters Pbi e, xa,xb and of the post- Keplerian parameter s. The mismodeling in them yields 



IPN 



{xa+xb) 



0.01845 deg yr" 



(^wipnI. < 2a;iPN (f ) = 0.01318 deg yr' 



IPN 



SX A 



(xa+xb) 



Se 



1 X 10 ^ deg yr 



2 X 10^6 degyr-\ 



(11) 



_ (^^ipnIp, < 3W1PN (^j = 0(10-8) degyr-i. 
Thus, the total uncertainty is 

(^wiPN < 0.03165 deg yr""^, 
which maps into a relative uncertainty of 

(5(jiPN 



1.8 X 10^^ 



WiPN 



(12) 



(13) 



'See also l lLvutikov and Thompsonl . [20051) . 

''in principle, also the IPN correction to the third Kepler law calculated by jPamour and Deruellelll98^ 
should be included, but it does not change the error estimate presented here. 
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As a consequence, we have the important result 

Aoj = Wmcas - wiPN = (-0.00463 ± 0.03233) deg yr'^ 



(14) 



Every attempt to measure or constrain effects predicted by known Newtonian and post- 
Newtonian physics (like, e.g., the action of the quadrupole mass moment or the gravitomag- 
netic field), or by modified models of gravity, for the periastron of the PSR J0737-3039A/B 
system must face with the bound of eq. ([T4l ). 

Should we decide to use both the post-Keplerian para meters related to the Shapiro delay 
( Damour and Deruelle . 19861 : Damour and Taylor . 1992 ) 



-2/3 



(15) 



for determining the sum of the masses, we would have, with eq. ([3]), 



WlPN 



(1 



2ir J 



9/4 



,19/4 



which yields 



oj^p^ = 17.25122 ± 2.11819 deg yr'^ 



(16) 



(17) 



The major source of uncertainty is r, with 2.06264 deg yr~^; the bias due to the other 
parameters is about the same as in the previous case. 

Let us, now, consider the second post-New tonian contribution to the periastron preces- 
sion dPamour and SchaefeA Il988l : IWexL Il995h 



W2PN 



3(G0JI)5/2 
c4a7/2(i _ e2)2 



13 f m\ + m| 



32 /mAmB\ 
^ 3 V A/2 ) 



(18) 



up to terms of order 0{e'^). For our system it amounts to 4 x 10~^ deg yr~^, so that it 
should be taken into account in Acj. However, it can be shown that the bias induced by 
the enws in 9Jl and a amounts to 4 x 10~® deg yr~^, thus affecting the gravitomagnetic 
precession at the percent level. 



2.2. Summary and Discussion 

O'Connelll (l2004h . aware of the presence of other non-gravitomagnetic contributions to the 



pulsar's periastron rate, proposed to try to measure the gravitomagnetic spin-orbit preces- 
sion of the orbital angular momentum (IBarker and O'Connelll. Il975al) (analogous to the 
Lense-Thining node precession in the limit of a test particle orbiting a massive body) 
which is not affected by larger gravitoelectric contributions. However, its magnitude is 
?a (10~^ deg yr~^) sin V', where il) is the angle between the orbital angular momentum 
and the pulsar's spin; thus, it would be negligible i n the PSR J0737-3 039A/B system be- 
cause of the near- alignment between such vectors (IStairs et al.,l2006l). in agreeme nt with 
the observed lack of profile variations ( Manchester et al. . 2005 : Kramer et al. . 2006h . 
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In regard to the measurement of the moment of inertia of the component A via the 
gravitomagnetic periastron precession, our analysis has pointed out that the bias due to the 
mismodelling in the IPN gravitoelectric contribution to periastron precession-expressed 
in terms of the phenomenologically measured parameters P^, e, xa,xb, s-is the most im- 
portant systematic error exceeding the expected gravitomagnetic rate, at present, by two 
orders of magnitude; the major sources of uncertainty in it are xb and s, which should be 
measured three orders of magnitude better than now to reach the 10% goal. The projected 
semimajor axis xa of A, if known one order of magnitude better than now, would induce 
a percent-level bias. Instead, expressing the IPN gravitoelectric periastron rate in terms of 
Ph,e,XA,XB,s,r would be definitely not competitive because the improvement required 
for r would amount to five orders of magnitude at least. We prefer not to speculate now 
about the size of the improvements in timing of the PSR J0737-3039A/B system which 
could be achieved in future. Since the timing data of B are required as well for .xb and 
in view of the fact that B appears as a strong radio source only for two intervals, each of 
about 10-min duration, while its pulsed emission is rather weak or eve n undetectable for 
most of the remainder of the orbit (ILyne et al.l . l2004l : lBurgay et al.l.l2005r) . the possibility of 
reaching in a neai" future the required accuracy to effectively constrain I a to 10% level or 
better should be considered with mor e skepticism than done so far. Another analysis on this 
topic has been recently performed by Kramer and Wex ( 20091 ') . 



3. Testing a Uniform Cosmological Constant and the DGP Grav- 
ity with the Pulsar J0737-3039A 

Since, at present, the only reason why the cosmological constanj^ A is be l ieved to be 
nonze ro relies upon the observed acceleration of the universe (|Riess et al.l.ll998l : IPerlmutter et al 
19991) . i.e. just the phenomenon for which A was introduced (again), it is important to find 
independent observational tests of the existence of such an exotic component of the space- 
time. 

Here we put on the test the hypothesis tha t A 7^ 0^ whe re A is the uniform cosmo- 
logical constant of the Schwarzschild-de Sitter dstuchlkl . E999I) (or Kottler dKottlej. [l918l) ) 
spacetime, by suitably using the latest determinations of the parameters (see Tabled]) of the 
double pulsar PSR J0737-3039A/B system. 

The approach followed here consists in deriving analytical expressions Oa for the ef- 
fects induced by A on some quantities for which empirical values Omcas determined from 
fitting the timing data exist. By taking into account the known classical and relativistic 
effects Oknown affecting such quantities, the discrepancy AO = Omeas — Oknown is con- 
structed and attributed to the action of A, which was not modelled in the pulsar data pro- 
cessing. Having some AO and Oa at hand, a suitable combination C, valid just for the case 
A 7^ 0, is constructed out of them in order to compare Cmcas to Ca- if the hypothesis A / 
is coiTcct, they must be equal within the enws. Here we will use the anomalistic period 
Pb and the periastron precession tuo for which purely phenomenological determinations ex- 
ist in such a way that our C is the ratio of Au to APy,; as we will see, this observable is 



see dCalder and Laha\il2008l) and references therein for an interesting fiistorical overview. 
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independent of A but, at the same time, it retains a functional dependence on the system's 
parameters peculiar to the A— i nduced force and of any other Hooke-like forces. 

This Section complements dloriol . 12008 al) in which a similar test was conducted in the 
Solar System by means of the latest determinations of t he sec i ilar pre cessions of the longi- 
tudes of the perihelia of several planets. The result by llorid (12008 al) was negative for the 
Schwarzschild-de Sitter spacetime with uniform A; as we will see, the same conclusion will 
be traced here in Section [3. 1.1.1 

A complementary approach to explain the cosmic acceleration without resorting to dark 
energy was follo wed by Dvali , Gaba dadze and Porrati (DGP) in their braneworld modified 
model of gravity ( Dvali et al. . 2000h . Am ong other thin gs, it predicts effects which could 
be tested on a local, astronomical scale. In (llorid. l2008al) a negative test in the Solar System 
was reported; as we will see in Section 1331 PSR J0737-3039A/B confirms such a negative 
outcome at a much more stringent level. 

The overview and the conclusions are in Section 1331 



3.1. The Impact of A on the Periastron and the Orbital Period 



The Schwarzschild-de Sitter metric induces an extra-acceleratioqj (|Rindlen. 1200 ih 



-Ac^r, 
3 



(19) 



where c is the speed of light; eq. (|47] ). in view of the extreme smallness of the assumed 



nonzero value cosmological constant (A « 10 m ^), can be trea ted perturbatively with 



the standard techniques of celestial mechanics. In (iKerr et al.L l2003l) the secular precession 
of the pericentre of a test particle around a central body of mass Tl was found to be 



where 



n 



cm 



(20) 



(21) 



is the Keplerian mean motion; a and e are the semimajor axis a nd the eccentrity, respec- 
tively, of the test particle's orbit. Concerning a binary system, in djetzer and Serend . 120061) 
it was shown that the equations for the relative motion are those of a test particle in a 
Schwarzschild-de Sitter space-time with a source mass equal to the total mass of the two- 
body system, i.e. Tl = xua + niB- Thus, eq. (l20l ) is valid in our case; a is the semi-major 
axis of the relative orbit. 



Following the approach by lJetzer and Serend (l200q) . we will now compute Pa, i.e. the 
contribution of A to the orbital period. One of the six Keplerian orbital elements in terms 
of which it is possible to parameterize the orbital motion in a binary system is the mean 
anomaly A4 defined as 7W = n{t — Tq), where n is the mean motion and Tq is the time 
of pericenter passage. The mean motion n = 2tt/P\^ is inversely proportional to the time 
elapsed between two consecutive crossings of the pericenter, i.e. the anomalistic period 



The present test is valid for all exotic Hooke-type forces of the form Cr i lCalder and L^iavlbOOSl) , with C 
arbitrary nonzero constant. 
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Pb- In Newtonian mechanics, for two point-like bodies, n reduces to the usual Keplerian 
expression n = 27r/Pb- In many binary systems, as in the double pulsar one, the period 
Pb is accurately determined in a phenomenological, model-independent way, so that, in 
principle, it accounts for all the dynamical features of the system, not only those coming 
from the Newtonian point-like terms, within the measurement precision. 

The Gauss equation for the variation of the mean anomaly, in the case of an entirely 
radial disturbing acceleration Ar like eq. (l47l ). is (IRoyi 119881) 



dM 2 ^ /r\ (l-e2) ^ 

— = n Ar - + ^ -Ar COS /, 

at na \a/ nae 



(22) 



where / is the true anomaly, reckoned from the pericenter. Using the eccentric anomaly E, 
defined as 

M = E - esinE, (23) 

turns out to be more convenient. The unperturbed Keplerian elUpse, on which the right- 
hand-side of eq. (|22] ) must be evaluated, is 



r = a (1 — e cos E) ; 



(24) 



by using eq. (|47] ) and 



eq. ((22)) becomes 



( ^ = l-ecosE, 



/( 
= T 



c os E—e 
-e cos E ' 



(25) 



dE 
dt 



n 



(1 — e cos E) 



Ac2 



2(1 - ecos-E) - 



(cos E 



(26) 



Since Ac^/3n^ « 10 from eq. (l26l ) it can be obtained 



^ ^'■^'^ (1 - ecos£;) r Ac^ 



2(1 - ecos£;)' 



(cos E — e) 



dE, 



(27) 



which integrated yields that 



with 



Pa 



Pb = Pb + Pa 

_ ttAc^ (7 + 3e^) 
~ 3^3 



(28) 



(29) 
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3.1.1. Combining the Periastron and the Orbital Period 



For the sake of convenience, from Section I2X] let us recall the general relativistic expres- 
sions of the post-Keplerian parameters r, s and ui; 



TqITIb, 



CJiPN 



1-e^ 



-2/3 



(30) 



-5/3 



By means of 



-{XA + xb) 



(31) 



and of the equations for r and s it is possible to express Py^ and cI^ipn in terms of Pb and of 
the phenomenologically determined Keplerian and post-Keplerian parameters XA,XB,r,s 
as 

1/2 



WiPN 



3/2 (XA-\^/^ g-9/4. 



(32) 



3sr 



XA{l-e') 



In such a way we can genuinely compare them to Pb and co because they do not contain 
quantities obtained from the third Kepler law and the general relativistic periasti'on preces- 
sion themselves; moreover, we have expressed the sum of the masses entering both Pb and 
"^iPN in terms of r and s, thus avoiding any possible reciprocal imprinting between the third 
Kepler law and the periastron rate. At this point it is possible to construct 



R 



Aw 
AP' 



with 



note that 




Ph 



i^lPN, 



R = R{P]„XA,XB,e;uj,r,s). 



(33) 



(34) 



(35) 



By attributing Aco and AP to the action of A, not modelled into the routines used to fit the 
PSR J0737-3039A/B timing data, it is possible to compare R to 



Ra 



UJA 



37r3^Pj,^3/2^9/2 



4^2(7 + 3g2)(^^ + 

xb^xJ 



(36) 



and see if eq. (1331) and eq. (1361) are equal within the errors. Note that eq. (1361 ) is independent 
of A and, by definition, is able to test the hypothesis that A 7^ 0. From Table [T]it turns out 



Pa = (3.4 ± 0.3) x 10"^ 



(37) 
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R\ is a well determined quantity, different from zero at about 1 1 sigma level. In regard to 
R we have 

( Au = -0.3 ±2.1 degyr-\ 

(38) 

AP = 59 ± 364 s, 



so that 



\R\ = (0.3 ±4) X 10"^^ s~ 



(39) 



R is compatible with zero in such a way that its range does not overlap with the one of ^a: 
indeed, the upper bound on R is three orders of magnitude smaller than the lower bound on 
Ra- Thus, we must conclude thaj^ 

R^Ra- (40) 

Concerning the released uncertainties in R and Ra, they must be considered as upper 
bounds since they have been conservatively computed by linearly adding the individual 
biased t erms due to 6Ph, Soj , 6e, 6xa, Sxb, 6r, 6s in order to account for the existing corre- 
lations (jKramer et al.l . l2006r) among them. 

The results of the present study confirm those obtained in the Solar System by taking 
the ratio of the estimated corrections to the standard Newtonia n/Eins t einian precessions of 
the longitude of the perihelia zu for different pairs of planets dloriol 12008 al) . It would be 
very interesting to devise analogous tests involvi ng other observables (lensing, t ime delay) 
affected by A as well recently computed in, e.g., ( Ruggiero . 2009; SerenoL 2008). 



3.2. The Dvali-Gabadadze-Porrati Braneworld Model 

The approach prev iously outlined for A can be followed also for the DGP braneworld model 



(IDvali et al.L 12000) which recently re ceived great attention from an observational point of 
view (lDvaUetal.l.l2003l : lloriol. l2008bh . 

The preliminary confrontations with data so far performed refer to the perihelia of the 
Solar System planets. Indeed, DGP gravity predicts a n extra-precession of the pericentre of 
a test particle ( Lue and StarkmanL 2003; loriol. 20053 ) 



wdgp = =Ft 



8 V ro 



1 



13 
32^ 



(41) 



where the signs =F are related to the two different cosmological branches of the model and 
ro is a free-parameter set to about 5 Gpc by Type lA Supernovse data, independent of the 
orbit's semimajor axis. The predicted precessions of about 10~^ arcsec cy~^ were found 
to be compatible with the estimated correc t ions to the usual apsidal precessions of planets 
considered one at a time separately ( loriol. 2008b|), but m arginally incompatible with the 

ratio of them for some pairs of inner planets (|lorio . 2007b ). 

The effects of DGP model on the orbital period is ( loriol 20061) 



-Pdgp 



11 



2\2 



a^l-e^) 

cm 



(42) 



'in principle, also the IPN correction to the third Kepler law jPamour and Deruellelll98^ should be in- 
cluded in AP, but it does not change the result. 
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From eq. (|4TI) and eq. (1421) it is possible to construct 

. '^DGP 



-Rdgp 



-Pdgp ' 



(43) 



which, expressed in terms of the phenomenologically determined parameters of PSR J0737 
3039A/B , becomes 

3 (1 _ 13 p^^^3/2^9/2 



R 



DGP 



22n{l-e'^){xx + XBfx^'^' 



(44) 



Putting the figures of Table eq. (O into eq. (|44l) and computing the uncertainty as done in 
the case of A yields 

i?DGP = (1.4 ± 0.1) X lO^'^ s"^ (45) 

As can be noted, the lower bound of -Rdgp is four orders of magnitude larger than the upper 
bound of R, so that we must conclude that, also in this case. 



R / -Rdgp. 



(46) 



The outcome by loriol ( 2007b ') is, thus, confirmed at a much more stringent level. 

An analysis of typ e la supernovae (SNe la) data disfavoring DGP model can be found 
in dfiento et al.ll2005h . 



3.3. Summary and Discussion 

In this Section we used the most recent determinations of the orbital parameters of the dou- 
ble pulsar binary system PSR J0737-3039A/B to perform local tests of two complementary 
approaches to the issue of the observed acceleration of the universe: the uniform cosmo- 
logical constant A in the framework of the known general relativistic laws of gravity and 
the multidimensional braneworld model by Dvali, Gabadadze and Porrati which, instead, 
resorts to a modification of the cun^ently known laws of gravity. Since, at present, there are 
no observational evidences for such theoretical schemes other than just the cosmological 
phenomenon for which they were introduced, it is important to put them on the test inde- 
pendently of the cosmological acceleration itself. It is worthwhile noting that the results for 
A hold also for any other Hooke-like additional force proportional to r. 

To this aim, we considered the phenomenologically determined the periastron pre- 
cession Co and the orbital period Pb of PSR J0737-3039A/B by contrasting them to the 
predicted IPN periastron rate o^ipn and the Keplerian period Pb- With such discrepan- 
cies we constructed the ratio R = Atj/AP by finding it compatible with zero: \R\ = 
(0.3 lb 4) X 10^^^ s^^. Then, we compared R to the predicted ratios for the effects of A and 
the DGP gravity-not modeled in the pulsar data processing-on the periastron rate and the 
orbital period by finding Pa = (3.4±0.3) x 10"^ s'^andPoGP = (1.4±0.1) x lO"'^ 
respectively. Thus, the outcome of such a local test is neatly negative, in agreement with 
other local tests recently performed in the Solai" System by taking the ratio of the non- 
Newtonian/Einsteinian rates of the perihelia for several pairs of planets. 



Some Applications of Binary Pulsars to Fundamental Physics 



13 



4. The IPN Secular Effects on the Mean Anomaly in Binary Pul- 
sar Systems 

According to the Einstein's General Theory of Relativity (GTR), the post-Newtonian gravi- 
toelectric two-bo dy acceleration of order 0(c~ ^) (IPN) is , in the post-Newtonian cent re of 
mass frame (see ( Damour and Deruelle , 1985h and, e.g., ( Portilla and Villareal , 2004 ) and 
references therein) 



a 



GE 



cm 



r + (r • v){4: - 2v)v 



(47) 

where r and v are the relative position and velocity vectors, respectively, G is the Newto- 
nian constant of gravitation, c is the speed of light, mi and m2 are the rest masses of the 
two bodies, 9Jt = mi + m2 and u = mi'm2/m? < 1. 

Let us recall from Section 13.1.1 that the orbital phase can be characterized by the mean 
anomaly Ai defined as 

M = n{t-To), (48) 



(49) 



where the unperturbed mean motion n is defined as 

. 27r 

In it Pf, is the anomalistic period, i.e. the time elapsed between two consecutive pericentre 
crossings, which is 2n / GTl for an unperturbed Keplerian ellipse, and Tq is the date of 
a chosen pericentre passage. 

The variation of the mean anomaly can be written, in general, a^ 



dM 



n 



27r 



(50) 



The second term of the right-hand side of eq. (1501 ) accounts for any possible variation of 
the anomalistic period. The third term, induced by any small perturbing acceleration with 
respect to the Newtonian monopole, whether relativistic or not, is the change of the time of 
the pericentre passage, which we will define as 



'dt~ Tb \ dt ) ■ 
It can be calculated with the aid of the Gaus^ perturbative equation ( Roy , 1988 ) 

2 



dt 



na 



duj dVl\ 
— + cosi—- 
dt dt J 



(51) 



(52) 



where i, are the inclination and the longitude of the ascending node, respectively, of the 
orbit. In order to obtain the secular effects, we must evaluate the right-hand-side of eq. 



'^See also lLin- Sei] l l2010l) . 

^F or a different approach based on a modified form of tfie Lagrange planetary equations see dCalura et aI.L 
Il997t) . 
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on the unperturbed Keplerian ellipse and, then, average the result over one orbital 
revolution. 

We will now consider eq. (|47] ) as perturbing acceleration. Let us start with the first term 
of the right-hand-side of eq. ( [52l ). By defining 



A = 


^(4 + 2.) 


B = 


-^(1 + 3^^), 


C = 


^(4-1), 



(53) 



it is possible to obtain from eq. ([47 



GE 



A 



+ B 



+ C 



(54) 



Now the term —2Arr / no?, with Ar given by eq. (l54l) . must be evaluated on the unperturbed 
Keplerian ellipse chai^acterized by 



a(l-e^) 
1+e cos / ' 

nae sin / 



^(l + e2 + 2ecos/) 
where / is the true anomaly, and averaged over one orbital period by means of 



(55) 



r'^df 



dt 



(56) 



Thus, 



na ^ a J \Pb J 



r\/l - 



[A + Brv"^ + Cr{rf] df. (57) 



In the expansion of r in eq. (l57l) the terms of order 0{e^) are retained. The final result is 



na a 



ncm 



Wl — e- 



.H{e-v), 



(58) 



with 



F ~ - 2(4 + 2i/) + (1 + 3z^) ( 2 + + — + 



(59) 
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The post-Newtonian gravitoelectric secular rate of pericentre is independent of v and is 
given by the well known formula 



Itt 



GE c^a{l-e^y 



(60) 



while there are no secular effects on the node. 

The final expression for the post-Newtonian secular rate of the mean anomaly is ob- 
tained by combining eq. (l58])-eq. (l60l ) and by considering that, for a two-body system, it is 
customarily to write 

nGTi _ /Pb}~^^^ 
~ \2^, 

It is 



{TqM) 



2/3 



(61) 



dt 



GE 



-9(£) {T^Mfl\l-e')-^l'F{e-u) 



with 



e2 e6 



9 



16 



32 



(62) 



(63) 



Note that eq. (1621 ) is negative because eq. (1631) is always positive; thus the crossing of the 
apsidal line occurs at a later time with respect to the Kepler-Newton case. 

Note that, for — )• 0, i.e. mi ^ m2, eq. (l62l ) does not vanish and, for small eccentrici- 
ties, it becomes 



di 

dt 



9nGm2 



GE 



(64) 



which could be used for planetary motion in the Solar System ( loriol . 2005b ). E.g., for 
Mercury it yields a secular effect of almost —130 aixsec cy~^. It is important to note that 
the validity of the present calculations has also been numerically checked by integrating 
over 200 years the Jet Propulsion Laboratory (JPL) equations of motion of all the plan- 
ets of the Sola r System with an d without the gravitoelectric 1/c^ terms in the dynamical 
force models (lEstabrookl Il97ll) in order to single out just the post-Newtonian gravito- 
electric effects. They fully agree with eq. (l64l ) (E.M. Standish, private communication, 
2004). Another analytical calculation of the post-Newtonian genera l relati vistic gravito- 
electric secular- rate of the mean anomaly was performed (Rubincam . 1977b in the frame- 
work of the Lagrangian p erturbative scheme for a central body of mass 9JT-test particle 



system. iRubincamI (|1977n starts from the space-time line element of the Schwarzschild 
metric written in terms of the Schwarzschild radial coordinate r . Instead, eq. (|47] ) and 
the equations of motion adopted in the practical planetary data reduction at, e.g. JPL, are 
written in terms of the standard isotropic radial coordinate r related to the Schwarzschild 
coordinate by r = r(l + (?5Jt/2c^r)^. As a consequence, the obtained exact expression 



dt 



(Rubincam) 



GE 



2,nGm 



(65) 



contrary to the pericentre case, agrees neither with eq. (1641 ) nor with the JPL numerical 
integrations yielding, e.g., a secular advance of +42 arcsec cy~^ for Mercury. For a better 
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und erstanding of such comparisons, let us note that both the numerical analysis by Standish 
and iRubincamI (119771) are based on the Ph = ca se; n gets canceled by construction in 
the Standish calculation, while in (lRubincaml . [l977l) the numbers are put just into eq. 
which is the focus of that work. 



4.1. Testing Gravitational Theories with Binary Pulsars 

In general, in the pulsar's timing data reduction proces^ five Keplerian orbital parameters 
and a certain number of post-Keplerian parameters are determined wi th great accuracy in 
a phe nomenological way, independently of any gravitational theory dWexl. |200l|; IStairsl. 
2003|). The Keplerian parameters are the projected semimajor axis x = as'mi/c, where i is 
the angle between the plane of the sky, which is normal to the line of sight and is assumed as 
reference plane, and the pulsai^'s orbital plane, the eccentricity e, the orbital period P5, the 
time of periastron passage Tq and the argument of periastron coq at the reference time Tq. 
The most commonly used post-Keplerian parameters are the periastron secular advance u, 
the combined time dilation and gravitational redshift due to the pulsar's orbit 7, the variation 
of the anomalistic period P^, the range r and the shape s of the Shapiro de l ay. Th ese post- 
Keplerian parameters are included in the timing models dWexL I2OOII : IStairsL |2003|) of the so 
called Roemer, Einstein and Shapiro Ar, Ae, As delay j^'l occurring in the binary pulsar 
systertJ^ 



Ar = 



Af 



;sina;[cos£' — e(l + 6r)] + x cos uj sin E 1 — 6^(1 + 5e 



■jsinE, 



(66) 



-2rln{l — ecosE — s[sina;(cos -E — e) + \/l — cos w sin E] } , 



where E is the eccentric anomaly defined as E — e sin E = A4. cos E and sin E appearing 
in eg. (l66l) c an be expressed in terms of Ai by means of the following elliptic expansions 
(IVintiLll998h 

( cosS = -I + ^Jl^f,£[J^{je)]cos{jM), 

(67) 

, sini? = l^-^i^sin(jA^), 
where Jj{y) are the Bessel functions defined as 



irJjiy) 



cos{j9 — y sin 



(68) 

can be accounted for in the 



The relativistic secular advance of the mean anomaly eq. 
pulsar timing modelling by means of eq. (|67] ). 

In a given theory of gravity, the post-Keplerian parameters can be written in terms of 
the mass of the pulsar mp and of the companion mc. In general, and nxc are unknown; 



For all general aspects of the binary pulsar systems see Jwe} ! I2OO ll ; fstairsl l2003h and references therein. 
"For the complete expression of the timing models including, e.g., also the delays occurring in the Solar 
System due to the solar gravity see ijWex, 2001; Stairs, 2003). 

'^The aberration parameters 5r and 5e are not, in general, separately measurable. 
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this means that the measurement of only one post-Keplerian parameter, say, the periastron 
advance, cannot be considered as a test of a given theory of gravity because one would not 
have a theoretically calculated value to be compared with the pheno menologically measured 
one. In GTR the previously quoted post-Keplerian parameters are (IDamour and Deruelle , 
1986h 



1 /3 

P - 192^^5/3 (P,Y^'^ (l+ie' + i^') mpr 
— 5 -^G \27rJ (l-e2)7/2 Ml 



/3 ) 



-2/3 



(69) 



It is important to note that the relativistic expression of Pb in eq. (1691 ). should not be con- 
of eq. ( [62l ). Indeed, it refers to the shrinking of the orbit due to grav- 



GE 



fused with ^ 

itational wave emission which vanishes in the limit v ^ 0, contrary to eq. ( [62] ) which 
expresses a different, independent phenomenon. The measurement of two post-Keplerian 
orbital parameters allows to determine mp and nXc, assumed the validity of a given theory of 
gravity ^^1. Such values can, then, be inserted in the analytical expressions of the remaining 
post-Keplerian parameters. If the so obtained values are equal to the measured ones, or the 
curves for the 2 + A'', with A'' > 1, measured post-Keplerian parameters in the mp — nXc 
plane all intersect in a well determined (mp,mc) point, the theory of gravity adopted is 
consistent. So, in order to use the pulsar binary systems as valuable tools for testing GTR 
the measurement of at least three post-Keplerian parameters is required. The number of 
post-Keplerian parameters which can effectively be determined depends on the character- 
istics of the particu lar binary system under consideration. For the pulsar-neutron star PSR 
B 19 13-1- 16 system (|Hulse and Tayloii . 1 197 50 the three post-Keplerian parameters lo,j and 
Ph were measured with great accuracy. For the pulsar-neutron star PSR B 1534-1-12 system 
(|Stairs et al.l . l2002|) the post-Keplerian parameters reliably measured are to, 7, r and s. For 
the pulsar-pulsar PSR J0737-3039 A-i-B system the same four post-Keplerian parameters 
as for PSR B 1534-1- 12 are available plus Pf, and a further constraint on nip/mc coming 
from the measurement of both the projected semimajor axes. On the contrary, in the pulsar- 
white dwarf binary systems, which are the majority of the binary systems with one pulsar 
and present almost circular orbits, it is often impossible to measure 00 and 7. Up to now, 
only r and g have been measured, with a certain accuracy, in the PSR B 1855-1-09 system 
(jKaspi et al.L 119941 ). so that it is impossible to use its data for testing the GTR as previously 
outlined. 



This would still not be a test of the GTR because the masses must be the same for all the theories of gravity. 



of course. 
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4.2. The Secular Decrease of the Mean Anomaly and the Binary Pulsars 



Let us investigate the magnitude of the mean anomaly precession in some systems including 

one or two radiopulsars. 

For PSR B 1913+16 we have (|Weisberg and Taylotl.l2005h ruj. = 1.4414, ruc = 1.3867,e 
0.6171338, Pb = 0.322997448930 d. Then, u = 0.2499064, F = 1.04459537192 and 
-10.422159 deg yr~^. For PSR J0737-3039 A we have u = 0.249721953643, 



GE 



0.946329857430. Thus, ^ 



GE 



-47.79 deg yr ^. This implies that the ratio of 



the post-Newtonian gravitoelectric secular rate of the mean anomaly to the mean motion 
amounts to ~ 10~^. Let us see if such post-Newtonian shift is detectable from quadratic 
fits of the orbital phases of the form M = qq + bot + CQt^. For PSR B1913-I-16 the 
quadratic advance due to the gravitational wave emission over thirty years amounts to 

(A 



-2.4184 X 10-12) 



-vr 



gw 



V 



{t - Tof = 0.5 deg, 



(70) 



with an uncertainty 6{AM) fixed to 0.0002 deg by 6Pb 
due to eq. ( l62l ) amounts to 



0.0009 X 10-^^. The linear shift 



GE 



(t - To) = -312.6647 deg 



(71) 



over the same time interval. The uncertainty ir0 n amounts to 1 x 10 ^ deg yr ^ due to 
(5Pb = 4 X 10-12 For PSR J0737-3039 A the gravitational wave emission over three years 
(Pb = —1.20 X 10-12) induces a quadratic shift of 0.008 deg, with an uncertainty 6{AM) 
fixed to 0.0005 deg by 6Pb = 0.08 x lO-^^. The hnear shift due to eq. ^ amounts 
to -143.3700 deg over the same time interval. The uncertainty in n amounts to 7 x 10-'' 
deg yr"! due to 6Pb = 2 x 10- ^'^ d. Thus, it should be possible to extract ^ from the 

GE 

measured coefficient bo; both the corrupting bias due to the uncertainties in the quadratic 
signature and the errors in n would be negligible. 



Measuring ^ 



GE 



as a further post-Keplerian parameter would be very useful in those 



scenarios in which some of the traditional post-Keplerian parameters are known with a 
modest precision or, for some reasons, cannot be considered entirely reliable^. E.g., in the 
double pulsar system PSR J0737-3039 A -i-B the parameters r and 7 are measured with a rel- 
atively low accuracy (ILyne et all 120041) . Moreov er, there are a lso pulsar binary systems in 
which only the periastron rate has been measured ( Kaspi . 1999h : in this case the knowledge 
of another post-Keplerian parameter would allow to determine the masses of the system, 
although it would not be possible to constraint alternative theories of gravity. 



''^The sum of the masses and the semimajor axis entering n are determined from timing data processing 
independently of M itself, e.g. from the periastron rate and the projected barycentric semimajor axis. 

'^The measured value of the derivative of the orbital period Pb is aliased by several external contributions 
which o ften limit the precision of the tests of competing theories of gravity based on this post-Keplerian pa- 
rameter ( IStairsll2003h . 
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4.3. Overview and Discussion 

In this Section we have analytically derived for a two-body system in eccentric orbits the 

secular variation ^ yielding the post-Newtonian general relativistic gravitoelectric part 

GE 

of the precession of the mea n anomaly not d ue to the variation of the orbital period. For a 



complementary analysis, see lLin-SenI (I2010f) . In the limit of small eccentricities and taking 
the mass of one of the two bodies negligible, our results have been compai^ed to the out- 
come of a numerical integration of the post-Newtonian general relativistic gravitoelectric 
equations of motion of the planets of the Solar System performed by JPL: the agreement 
between the analytical and numerical calculation is complete. Subsequently, we have in- 
vestigated the possibility of applying the obtained results to the binary systems in which 
one pulsar- is present. In particular, it has been shown that the variation of the orbital period 

Pi) due to gravitational wave emission and the effect derived by us are different ones. 

gw 

Indeed, the post-Newtonian gravitoelectric precession of the mean anomaly, which is al- 
ways negative, is related to the secular increase of the time of pericentre passage and occurs 
even if the orbital period does not change in time. A quadratic fit of the orbital phase of the 
pulsar would allow to measure £ because the biases due to the errors in the quadratic 

GE 

shift due to Pi, and in the Unear shift due of the mean motion n are smaller. The use of 
^ as a further post-Keplerian parameter would allow to improve and enhance the tests 

GE 

of post-Newtonian gravity especially for those systems in which only few post-Keplerian 
parameters can be reliably measured. 
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